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^ _  INTRODUCTION 

An  historical  and  rather  complete  statement  of  the  atmospheric  ozone 
problem  has  been  given  by  Ootz  (1951)  and  Craig  (1950)  and  hence  such  a  state¬ 
ment  in  this  writing  would  serve  no  purpose.  It  is  well  to  tabulate  the 
observations  and  factors,  however,  which,  prior  to  this  writing,  were  considered 
important  to  the  complete  ozone  problem. 

1.  There  is  a  seasonal  variation  (spring  maximum  and  fall  minimum) 
of  the  total  ozone  in  a  vertical  column.  Apparently  there  is  also  a  late 
summer  maximum  observed  at  all  latitudes. 

2.  There  is  a  latitudinal  variation  of  the  total  ozone  in  a  vertical 
column  with  a  maximum  at  45°  -  6o°  N  latitude  and  a  minimum  in  the  equatorial 
regions  throughout  the  year.  A  southerly  latitude  maximum  is  not  clearly 
defined,  however  observations  in  the  southern  hemisphere  are  sparse. 

3.  Large  day  to  day  variations  exist  which  are  greater  than  the 
yearly  variation  of  the  monthly  averages.  It  has  been  suggested  that  latitu¬ 
dinal  adveotion,  vertical  motion,  solar  production  end  dissociation  alone  or 
in  combination  are  responsible  for  these  large  day  to  day  variations. 

4.  The  vertical  distribution  of  ozone  has  not  been  explained  satis¬ 
factorily  although  the  measurements  and  the  photochemical  theory  show  a  middle 
stratosphere  maximum. 

5.  The  Umkehr  measurements  show  a  relationship  between  the  vertical 
distribution  and  the  total  ozone  in  a  vertical  column.  The  Umkehr  technique 
seems  to  suggest  that  the  major  variation  in  the  total  ozone  occurs  between  10 
and  20  km.  It  also  appears  that  there  is  a  positive  correlation  between  the 
potential  temperature  at  12,  15  and  18  km  and  the  total  ozone,  while  there  is 
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a  negative  correlation  between  the  tropbpause  height  and  the  total  ozotie. 

6.  A  correlation  between  total  ozone  and  solar  visil:|fility  (sun 

ii 

spot  number)  is  not  well  established  even  with  regard  to  sign. 

It  should  be  noted  that  in  the  above  ho  mention  is  made  of  the  sink  at 
ground  levels.  All  emphasis  has  been  given  to  variations  in  the  production 
rate  of  ozone  and  transport  phenomena  within  the  atmosphere.  This  is  unfortu¬ 
nate  since  the  variation  in  the  sink  strength  at  ground  level  may  well  doninate 
the  variation  in  production  at  high  levels.  Moreover,  in  the  above  tabulation, 
sources  of  ozone  other  than  the  production  from  solar  ultraviolet  light  have 
been  omitted.  It  will  be  shown  later  that  the  production  of  ozone  in  the 
troposphere  by  electrical  discharge  cannot  be  ruled  out  as  being  insignificant. 

In  this  writing  it  will  bo  the  main  concern  of  the  author  to  discuss  the 
vertical  distribution  of  ozone,  the  sink  at  ground  levels  which  removes  ozone 
from  the  atmosphere,  vertical  transport,  ozone  production  via  electrical  dis¬ 
charges  and  measurements  of  the  foregoing,  as  well  as  a  brief  explanation  of 
the  development  of  the  chemiluminescent  detector  applied  to  the  measurements. 
Lastly,  a  report  on  some  measurements  of  the  small  ion  number  within  the  free 
atmosphere  from  ground  level  to  balloon  altitudes  and  the  results  of  such 
measurements  will  be  given. 

The  discovery  of  the  chemiluminescent  reaction  between  ozone  and  solutions 
of  different  fluorescent  and  non- fluorescent  d3res  is  not  new  (Kearney,  1924j 
Biswas  and  Dhar,  1928) .  The  chemiluminescent  hydrazldes  in  alkaline  solution, 
particularly  lumlnol,  give  high  chemiluminescent  efficiencies  (Harvey,  1929} 
Briner,  1940;  Wilhelmsen,  Lurary  and  Eyring,  1955).  Kautsky  and  Zocher  (1921) 
were  original  in  obtaining  a  "dry"  chemiluminescent  reaction  between  ozone  and 
Rhodamine  B  adsorbed  on  slloxene.  Bernanose  and  Rene  (1959),  used  luminol 
for  ozone  determinations  by  adsorbing  it  on  chromatographic  paper  from  basic 
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solutions  and  exposing  the  paper  to  gaseous  ozone.  The  light  output  vas 
measured  with  a  phctoniultiplier.  It  was  not  until  196l»  however,  that 
V.  Rogener  brought  this  technique  to  the  foreground  and  applied  it  to  the 
measurement  of  atmospheric  ozone.  At  this  point  Kroening  and  Ney  (1962) 
began  the  development  of  a  chemiluminescent  ozone  detector  Independently 
of  the  work  of  Regener. 


I 
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THE  DETECTOR 

A  majority  of  the  author's  work  with  chemilumlnesoent  ozone  detectors 
Involved  the  development  of  the  lurainol  silica  gel  detector.  Luminol  (5  amlno- 
2»  3-dihydrophthalazlne-l,  4-dione)  dissolved  in  alkaline  solution  is  adsorbed 
on  silica  gel.  The  resulting  complex  has  a  very  high  efficiency  in  the  chemi¬ 
luminescent  decomposition  of  ozone.  The  efficiency  of  the  reaction  increased 
with  increasing  pH  to  the  point  where  the  silica  gel  decomposed  and  reverted 
to  sodium  silicate  (water  glass).  The  base  used  was  sodium  hydroxide.  The 
reaction,  however,  was  drastically  sensitive  to  water  vapor.  Much  work  was 
done  to  control  this  reaction.  Attempts  to  saturate,  fix,  or  render  the  silica 
surface  hydrophobic  failed  in  the  sense  that  the  detector  was  never  really 
stable  enough  for  use  In  very  high  altitude  measurements.  The  attractive  high 
efficiency  of  the  reaction,  however,  leads  one  to  investigate  every  conceivable 
method  that  might  stabilize  the  reaction.  After  this  exhaustive  effort  failed, 
work  was  initiated  on  the  Rhodamlne  B  silica  gel  complex.  Rhodamine  B  was 
chosen  rather  than  other  fluorescent  or  nonfluorescent  dyes  because  the  litera¬ 
ture  Indicates  it  as  being  second  in  line  of  efficiency.  The  remainder  of  the 
detector  development  is  restricted  to  Rhodamlne  B,  not  because  the  ideas  to  be 
presented  would  not  apply  to  the  use  of  other  dyes,  but  rather  that  the  author 
worked  only  with  Rhodamlne  B,  Although  the  precise  nature  of  the  reaction 
between  ozone  and  Rhodamlne  B  adsorbed  on  silica  gel  is  not  clearly  understood, 
certain  features  of  the  reaction  are  evident  from  experiment.  The  reaction  is 
heterogeneous  since  it  takes  place  on  the  silica  surface.  It  is  not  clear 
whether  the  reaction  is  unlmolecular  or  blmolecular,  but  this  really  has  secon¬ 
dary  significance  for  its  use  in  ozone  measurements.  It  is  clear  from  experi¬ 
ments,  however,  that  the  reaction  is  linear  with  respect  to  the  gaseous  ozone 
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ooncentratldn  over  a  range  of  at  least  1000,  The  tnathematlcal 'formulation 

:! 

of  a  first  order  reaction  is  ;! 


(1)  to  ,  .wo  - 

vdiere  c  is  the  concentration  of  the  reactant  of  interest  (ozone).  The  oi?der 
of  a  reaction  is  determined  by  observation.  The  above  reaction  rate  is  easily 
understood  for  homogeneous  gas  reaction.  However,  its  application  to  a  hetero¬ 
geneous  surface  reaction  requires  a  further  linear  relationship  between  the  gas 
phase  concentration  and  the  surface  concentration  of  the  reactant.  This  linear 
relationship  supposedly  applies  in  the  case  of  the  ozone  Rhodamine  B-silica  gel 
complex  used  in  the  measurements  reported  in  this  work.  The  activation  energy 
required  for  the  chemiluminescent  decomposition  of  ozone  on  the  above  complex 
is  not  known.  Suffice  it  to  say  in  this  respect,  though,  that  under  the  temper¬ 
ature  conditions  in  which  the  detector  is  used  (0°C  -  30®C),  the  reaction  is 
not  significantly  temperature  sensitive.  The  silica  surface  acts  as  a  positive 
catalyst  in  the  chemlluminesenet  reaction  of  ozone  with  the  fluorescent  dye 


Rhodamine  B.  The  silica  gel  provides  a  polar  surface  upon  which  the  ozone  and 


Rhodamine  are  adsorbed  and  on  which  they  subsequently  react.  This  is  not  unlike 


the  manner  in  which  alcohols  and  basic  solutions  provide  surroundings  for  the 
same  reaction. 

Before  proceeding  with  the  chemistry  of  the  detector,  it  is  perhaps  well 
that  the  construction  and  use  of  the  detector  is  outlined.  Rhodamine  B  is  dis¬ 
solved  in  either  methyl  or  ethyl  alcohol.  Alcohols  of  higher  molecular  weight, 
except  for  propyl  alcohol,  polyhydroxyl  alcohols  such  as  glycerine  and  ethylene 
glycol,  are  unsuitable  for  chemilumlnescenoe  despite  the  fact  that  they  form 
beautiful  red-orange  fluorescent  solutions  with  Rhodamine  B.  This  may  possibly 
be  caused  by  steric  hindrance.  In  contrast  to  alcohols  and  alkaline  solutions, 
one  may  also  use  glacial  acetic  acid  as  a  solvent  in  which,  supposedly,  the 
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Hhodamine  B  becomes  assoolated  with  the  carboxyl  group,  When  Rhodamlne  B  in 
proper  solution  is  contacted  with  ozone*  light  Is  emitted  which  Is  not  greatly 
different  from  the  fluorescent  spectrum.  Nor  Is  the  Intensity  of  the  light  In 
the  case  of  methyl  and  ethyl  alcohol  or  acetic  acid  very  different.  However* 
for  polyhydroxyl  alcohols  or  alcohols  of  high  molecular  weight*  the  Intensity 
Is  considerably  diminished  or  undetectable.  It  Is  also  Interesting  that  the 
glow  Is  considerably  diminished  In  aqueous  solution. 

In  preparation  of  a  solution,  .05  gm  of  Rhodamlne  B  Is  dissolved  In  50  co 
of  solvent.  In  preparation  of  the  gel*  2  to  4  cc  of  solution  are  added  to 
2  ^  of  finely  divided  silica  gel  (200  mesh),  sold  commercially  as  acid 

t)y  Malllnckrodt  Chemical  Works.  The  type  of  silica  gel  used  Is  not 
Important  as  long  as  It  la  finely  ground  (high  specific  surface  area)  and  not 
contaminated.  The  gel  is  then  dried  in  a  vacuum  at  room  temperature.  To 
complete  construction  of  the  detector*  the  gel  Is  glued  onto  an  aluminum  disc 
(2"  dla  1/16"  thick).  The  disc  Is  then  placed  I/I6"  from  an  end  window  RCA 
6199  photomultiplier  with  proper  light  tight  fittings.  The  gas  to  be  sampled 
Is  admitted  to  the  small  space  between  dlso  and  photomultiplier  via  light  tight 
brass  tubing,  coated  on  the  Inside  with  flat  black  Krylon  paint.  After  evapor¬ 
ation  of  the  Krylon  propellant*  a  black  acrylic  resin  remains  on  the  tubing 
which  does  not  decompose  ozone  appreciably.  The  spacing  between  the  disc  and 
the  photomultiplier  Is  Important  In  order  to  obtain  the  maximum  signal  for  a 
given  ozone  concentration  and  yet  be  consistent  with  the  pumping  apparatus. 

The  spacing  of  I/I6"  was  found  to  optimum.  The  gas  to  be  sampled  is  aspirated 
past  the  disc  at  a  rate  of  approximately  500  cc  per  min.  At  very  high  flow 
rates*  the  signal  Is  nearly  Independent  of  flow  rate  where  evidently  an 
efficiency  factor  comes  Into  play.  At  a  flow  rate  of  ^00  cc  per  min,*  practi¬ 
cally  all  of  the  ozone  Is  removed  from  the  stream.  Choosing  a  pumping  speed 
such  that  all  the  ozone  is  removed  prevents  the  possibility  of  the  reaction  rate 
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being  controlled  by  diffusl,6n  to  the  surface.  A  gear  pump  driven  by  either  a 
oentrifugally  governed  d.e.  motor  or  a  synchronous  a.c.  motor  provides  the 
constant  fl^w  rate  required  for  the  measureioents.  The  pump  is  placed  on  the 
downstream  side  of  the  detector. 

The  photomultiplier  is  generally  operated  near  1200  volts  in  order  to 
provide  the  maximum  signal  as  veil  as  optimum  signal  to  noise  ratio.  A  22 
megohm  bleeder  network  provides  voltage  to  the  dynodes  atul  consists  of  2.2 
megohm  l/4  watt  resistors  between  dsrnodesi  except  between  the  cathode  and  first 
dynode  where  an  optimum  signal  to  noise  ratio  warrants  usirjg  3*9  megohms.  The 
photomultiplier  is  also  magnetically  and  electrostatically  shielded.  All  high 
voltage  sections  and  the  signal  leads  of  the  photomultiplier  are  potted  in 
silastic  rubber  to  prevent  high  voltage  breakdown  and  leakage  between  cathode 
and  anode. 

We  will  now  return  for  a  moment  to  the  chemistry  of  the  detector.  Two 
methods  were  used  to  verify  that  the  reaction  between  Rhodamlne  B  adsorbed  on 
silica  gel  and  ozone  in  gaseous  form  is  linear.  The  first  method  Involves  the 
use  of  an  ozone  generator  with  a  variable  but  calibrated  output.  The  detector 
is  subjected  to  the  variable  ozone  concentration  and  the  linearity  of  response 
is  quickly  checked.  The  second  technique  involves  placing  the  detector  in  an 
atmosphere  with  a  steady  but  arbitrary  ozone  concentration  and  then  varying  the 
total  pressure,  which  can  be  read  easily  with  a  manometer.  The  partial  pressure 
of  ozone  then  varies  in  proportion  to  the  total  pressure.  Such  a  test  has  been 
carried  out  in  a  bell  jar  which  was  connected  to  a  vacuum  pump.  An  initial 
ozone  level  is  established  and  then  decreased  in  proportion  to  the  total  pires- 
sure  by  the  vacuum  pump.  In  this  manner,  the  detector  has  been  checked  for 
linearity  over  a  range  of  partial  pressures  of  ozone  extending  from  about  1  m> 
mb  up  to  about  1000  ^  mb.  Presumably  the  detector  is  linear  over  a  much  greater 


range,  l?ut  the  need  to  verify  tbis  has  not  arisen.  Linearity  checks  using  thf 
first  technique  have  also  been  harried  out  and  yield  the  same  results.  There 
is,  however,  one  problem  with  the  detector  that  has  not  been  eliminated  comi^le- 
tely,  that  is,  its  sensitivity  to  water  vapor.  Yet,  this  is  not  surprising  when 
one  reflects  upon  the  universal  role  that  water  plays  in  chemical  reaotiofis .  It 
is,  so  to  speak,  the  universal  solvent  without  which  the  majority  of  chemical 
reactions  would  be  impeded,  if  not  prevented.  The  sensitivity  of  the  silica 
gel-Rhodamlne  complex  apparently  changes  by  25?^  to  100^  when  the  partial  pres¬ 
sure  of  water  is  changed  from  its  vapor  pressure  at  -72 up  to  a  concentra¬ 
tion  corresponding  to  70^  relative  humidity  at  25*C.  The  cause  of  the  variation 
in  sensitivity  from  one  detector  to  the  next  for  identical  variations  in  water 
vapor  has  yet  to  be  elucidated.  Much  effort  without  complete  success  has  been 
given  to  this  problem.  A  considerable  improvement,  on  the  average,  has  been 
achieved  by  using  a  substance  called  estersll  in  place  of  ordinary  silica  gel. 
Silica  gel  is  very  hydrophilic  and  is  used  as  a  drying  agent,  Estersll,  on  the 
other  hand,  is  a  silica  gel  which  has  been  esterified  (Her,  1955)*  The  polar 
end  of  the  molecule  is  permanently  attached  to  the  silica  surface,  leaving  the 
non-polar  and  extending  away  from  the  surface.  The  resultant  surface  is  so 
hydrophobic  that  immersion  in  water  does  not  wet  the  surface.  Nontheless, 
there  exists  on  the  average  a  variation  in  the  sensitivity  of  the  detector  of 
about  1/2  of  the  variation  given  for  the  ordinary  silica  gel.  The  use  of  a 
dessicant  ahead  of  the  detector  to  remove  the  water  vapor  immediately  suggests 
itself.  This  can  be  done  with  good  results  if  small  amounts  of  phosphorus 
pentoxide  dessicant.  The  disadvantage  of  other  des¬ 

sicants  is  that  they  decompose  the  incoming  ozone  as  well  as  remove  the  water 
vapor.  The  length  of  time  that  the  ^2©^  will  continue  to  function  as  a  dessi¬ 
cant  depends  upon  the  amount  of  ^2^5  present  and  the  total  water  vapor  influx. 
This  can  be  allowed  for.  The  other  alternative  is  to  select  the  detectors  with 
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the  least  variation.  The  ^.atter  alternative  has  been  .chosen  by  the  author. 
The  measurements  of  the  vertical  oaone  distribution  will  not  be  seriously 
affected  because  of  the  low  water  vapor  content  of  the  atmosphere,  except 
at  very  lo^V levels,  and  this  can  be  allowed  for  if  the  case  warrants.  The 

'I 

water  vapoV  problem  should  nontheless  be  investigated  farther  in  the  hope  that 
it  can  be  ^Alaiinated  eventually. 
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THE  OZONE  GENERATOR 

All  of  the  relative  measurements  which  can  be  made  with  this  ehemilumine- 
soent  detector  cannot  be  put  on  an  absolute  scale  until  one  has  an  absolute 
standai^  reference  source  of  osone.  To  date,  this  need  has  not  been  fulfilled. 

A  technique  using  long  lived  radioactive  material  5.s  suggested  from  the  measure¬ 
ments  of  Kroening  and  Nay  (1962),  but,  at  the  same  tine,  these  measurements 
indicate  that  the  radiation  levels  necessary  for  sufficient  ozone  production 
lie  well  in  the  lethal  range.  As  an  alternative,  a  secondary  ozone  standard 
has  been  built,  using  a  mercury  vapor  lamp  with  a  quartz  envelope.  This  secon¬ 
dary  standard  requires  periodic  calibration  by  means  of  an  absolute  chemical 
method  to  be  outlined  later.  The  disadvantage  of  the  secondary  standard  is 
that  only  its  average  output  over  a  period  of  time,  depending  upon  the  sensiti¬ 
vity  of  the  absolute  chemical  technique,  can  be  determined.  Short  term  stabi¬ 
lity  is  then  furnished  by  the  chemiluminescent  detector  to  calibrate  the  genera¬ 
tor  at  lower  levels  comparable  to  those  existing  in  the  atmosphere. 

A  description  of  the  secondary  standard  ozone  generator  is  now  presented. 
The  most  convenient  and  most  stable  generator  constructed  by  the  author  Uses  a 
General  Electric  4  watt  mercury  vapor  lamp  having  a  1  l/2"  diameter  quartz 
envelope.  These  bulbs  are  readily  available  and  are  used  mainly  for  odor  con¬ 
trol.  The  quartz  envelope  bulb  is  epoxied  to  a  polythylene  hat  such  that  only 
one  half  of  the  bulb  is  covered  by  the  hat  with  a  vol\ime  of  approximately  3  cc 
enclosed  by  the  hat  and  bulb  surfaces.  Into  this  volume,  one  Inserts  an  inlet 
and  outlet  through  which  the  air  to  be  ozonized  is  passed.  In  general,  care 
should  be  taken  with  all  surfaces  which  are  to  come  in  contact  with  the  ozone 
in  order  to  reduce  undesirable  decomposition  on  any  contaminated  surface. 


An  important  feature  of  the  generator,  however,  is  its  tendency  to  oleahse 
itself.  In  a  manner  similar  to  which  electronic  tubes  in  critical  circuits 
are  aged  for  stabilization,  so,  in  the  construction  of  the  generator,  the  bulb 
and  generating  chamber  are  aged  and  cleansed  merely  by  allowing  the  generator  to 
operate  for  a  considerable  length  of  time .  Two  days  of  continued  operation 
generally  suffice  to  attain  stability.  Yet,  prior  to  each  use  of  the  generator, 
a  short  warm-up  period  is  required.  The  output  of  the  generator  is  controlled 
the  power  supplied  to  the  bulb.  This  is  monitored  by  an  ammeter  and  -volt¬ 
meter.  An  autotransformer  provides  the  variation  of  bulb  current  and  voltage, 
and,  consequently,  the  output  of  the  generator.  Ballast  is  required  in  series 
with  the  bulb.  The  output  level  of  the  above  described  generator  generally  can 
be  varied  from  5  M'Xib  to  500  iinib  if  the  flow  rate  is  approximately  500  cc  per 
min  end  the  carrier  gas  is  air  at  NTP. 
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calibration  OF  THE  OZONE  GENERATOR 

The  oalibratlon  of  the  generator  described  in  the  previous  chapter  is 
accomplished  by  the  use  of  the  absolute  chemical  technique  described  Wadelin 
(1957) .  The  ozone  is  passed  through  a  solution  of  potassium  iodide  and  sodium 
thiosulfate.  The  iodine  is  displaced  as  follows: 

(2)  0^  +  2KI  +  HgO  — *>1^  Og 

The  iodine  thus  displaced  reacts  immediately  with  the  sodium  thiosulfate  accor¬ 
ding  to 

(3)  Ig  +  ZNagSgO^— ^HagSgOg  +  2NaI 

The  main  reason  for  this  step  is  to  prevent  the  volatization  and  subsequent  loss 
of  the  free  iodine  displaced  by  the  ozone.  The  above  solution  is  buffered  at 
pH  7  or  higher  in  order  to  assure  the  formation  of  no  more  than  one  mole  of 
iodine  per  mole  of  ozone.  At  low  pH,  one  would  have  the  side  reaction 

(4)  40.^  +  lOKl  +  lOK'*' - ►5I2  +  HgOg  +  4H2O  +  3O2  +  10K+ 

competing.  This  is  eliminated  by  the  buffer  solution.  In  equations  (2)  and 
(3),  the  reaction  is  complete,  i.e.,  proceeds  far  to  the  right. 

Since  potassium  iodate  (KIO^)  is  sufficiently  pure  to  be  used  as  a  stand¬ 
ard  and  is  stable  in  solution,  it  is  now  used  to  back  titrate  the  amount  of 
sodium  thiosulfate  remaining,  the  reaction  being 

(5)  10^-  +  51-  +  6H*— >312  +  3H2O 

Prior  to  the  back  titration  with  the  potassium  iodate,  a  quantity  of  sulfuric 
acid  is  added  so  that  the  above  reaction  will  proceed  far  to  the  right.  A 
blank  run  (no  ozone)  is  then  made  on  an  identical  sample  to  provide  the  blank 
sodium  thiosulfate  concentration.  A  subtraction  of  the  run  made  with  the 
generator  from  the  blank  then  yields  the  amount  of  sodium  thiosulfate  reacted 


with  the  iodine  which  was  displaced  by  the  ozone. 

The  end  point  in  the  above  titrations  is  determined  amperometrioally  using 
a  gold  and  silver  amperometrilc  probe  (Qergen,  1959) .  The  current  from  the  probe 
is  sent  through  a  variable  resistance.  The  voltage  across  this  resistance  is 
then  measured  with  a  0  -  10  mv  Leeds  and  Northrup  potentiometer.  The  variable 
resistance  (0  -  50K)  provides  a  wide  range  of  sensitivity.  In  this  way,  the 
end  point  is  plotted  automatically  by  the  recorder.  The  ozone  concentration 
is  computed  from  the  following  equations 

/XN  „  (A-B)  X  N  X  6  X  22.42  3C  760  X  T  X  10^ 


t^ere 


2xFxtxpx  273 
A  =  ml  of  KIO^  required  for  the  blank 
B  =>  ml  of  KIO^  required  for  the  ozone  sample 


N  =  normality  of  the  KIO^ 

T  =  ambient  air  temperature 
P  =  flow  rate  in  liters/min 
t  =  sampling  time  in  minutes 
p  ®  ambient  pressure  in  mm  Hg 

Z  =  ozone  concentration  in  parts  per  hundred  million  (pphm) 
of  air  by  volume.  The  procedure  used  in  the  above  analysis  is  given  in  an 
appendix. 
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7  ,  . 

OPERATHKl  PRpCEDURE 

TVie  ozone  measureTnents  being  reportesi  in  this  writing  are  basically  two- 
foldi  1.  continuous  measurements  made  1  meter  above  the  ground  in  Kinne- 

li 

apoliaj  2.  vertical  profiles  obtained  via  balloon  flights  launched  near  Mlnne- 

if  I . 

spoils.  The  ground  level  measurements  were  made  over  a  period  extending  from 
12/1/61  to  4/1/62.  Vertical  profiles  of  ozone  concentration  were  made  spora¬ 
dically  on  the  dates  indicated  in  the  graphs.  The  ground  level  monitor  consi¬ 
sted  of  a  Rhodamine  B-estersll  complex  and  photomultiplier  as  previously  out¬ 
lined  and  a  500  -  I600  V  high  voltage  power  supply  made  by  the  John  Fluke 
Manufacturing  Co.,  Seattle,  Washington.  A  Leeds  and  Morthrup  0  -  10  mv  poten¬ 
tiometer  with  a  chart  paper  drive  of  1  1/2**  per  hour  served  as  the  recorder. 

The  overall  sensitivity  of  the  system  was  checked  on  a  nearly  day  to  day  basis 
by  an  ozone  generator  of  the  type  previously  described.  Adjustments  were  made 
dally  to  maintain  the  sensitivity  of  the  system  such  that  a  full  scale  chart 
reading  corresponded  to  an  ozone  concentration  of  5  parts  per  hundred  million 
(5pphm)  of  air  by  volume,.  This  was  accomplished  by  adjusting  the  high  voltage 
to  the  photomultiplier.  Throughout  this  4-months  period,  the  sensitivity  of  the 
detector  decreased  by  about  a  factor  of  5  because  of  the  oxidation  of  the 
Rhodamine  and  contamination  of  the  detector  by  air  pollutants. 

The  detectors  used  in  the  high  altitude  soundings  were  of  the  same  type  as 
that  used  in  the  ground  monitor.  The  type  of  high  voltage  power  supply  used  in 
the  soundings  is  indicated  schematically  in  figure  1.  These  supplies  satis¬ 
fied  strict  stability  specifications.  The  gear  pump  is  driven  by  a  centri- 
fugally  governed  d.c.  motor,  providing  a  constant  volumetric  flow  rate  inde¬ 
pendent  of  the  pressure.  The  complete  flight  detector  weighs  only  1  1/2  lbs. 

The  detector  output  is  designed  to  work  dlxeetly  into  the  standard  I68O  me 
Weather  Bureau  radiosonde  system. 


PHOTOCHEKISTP.Y  OF  ATMOSPHERIC  OZONE 

Since  the  seasonal  Tariation  In  the  total  atmospheric  ozone  in  a  vertical 
column  is  180®  out  of  phase  with  the  sun,  little  hope  remains  for  the  photo¬ 
chemical  theory  alone  in  explaining  this  observation.  Transport  phenomena  and 
other  sources  of  ozone  are  undoubtedly  necessary  in  the  ocsnplete  explanation. 
Before  considering  other  processes,  it  is  convenient  to  consider  the  photo¬ 
chemical  theory.  The  photochemical  theory  originated  with  Chapman  (1930),  vras 
followed  up  by  Wulf  and  Deming  (1936),  and,  finally,  Craig  (1950)  has  given  much 
consideration  to  the  problem.  A  brief  sketch  of  the  theory  is  presented  below. 

The  dissociation  of  the  oxygen  molecule  is  taken  as  the  primary  step  in 
the  eventual  formation  cf  atmospheric  ozone.  This  proceeds  according  to  the 
following  scheme  discussed  originally  by  Chapman. 

a)  Og  +  hU-»0  +  0  l-sZiWO  A® 

b)  O  +  O  +  X-^Og  +  X 

c)  0  +  Og  +  X  — >  0^  +  X 

d)  0  +  0--e20^ 

3  2 

e)  0^  +  hP  — *  Og  +  0  X  <  11,000  A® 

The  Importance  of  three  body  collisions  relative  to  the  two  body  colli¬ 
sions  decreases  as  one  goes  to  lower  pressures.  That  is,  the  processes 

f)  0  +  0  — ^Og  +  hP 

g)  0  +  Og — »0^  +  hi> 

will  dominate  b)  and  c)  high  in  the  atmosphere. 

Let  n^,  Hg,  ny  and  n  be  the  0^,  0^,  0^  and  air  concentrations  respec¬ 
tively.  Denote  the  recombination  coefficient  for  processes  b),  c)  and  d)  by 
k^,  kj^,  and  k^^^  respectively.  Then  we  may  write  the  rate  equations  for 


•I 


i ' 
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and  n^  ag  foTioyts: 

h 


(7) 

(8) 


lit  ~  “*■  '^3"3  "  "  ^13’^”3 

*^''12W~\3V3 


& 

St 


'^3”3 


where  Jg  and  represent  quantities  proportional  to  the  solar  radiation 
available  for  processes  a)  and  e)  respectively. 

1  In  the  main  part  of  the  ozone  layer  process  d)  is  much  smaller  than 
process  c).  This  is  visualized  by  the  ratio  of  the  appropriate  terms  as 
follows  I 


(9) 


”3 


ki3nin3 


Bucken  and  Patat  (1936)  have  measured  the  ratio  which  is  repro¬ 
duced  in  figure  2,  Choosing  T  =  -80 ®C  and  n  —  3  X  10^®  equation  (9)  is  veri¬ 
fied,  Furthermore!  Chapman  (1955)  has  given  the  values  of  and  k^j^  as 

-  5  X  10"^^  cm^/s 


(10) 


^sec 


k^^  5  X  10"^  cm^/sec 


It  should  be  kept  in  mind  that  they  are  only  approximate  values  as  yet  to  be 
verified  by  experiment.  Using  these  values  one  likewise  finds  that  term  b)  is 
negligible  in  the  main  ozone  layer.  Setting  SJgUg  =  SQg  and  J^n-j  =  8et 

from  equation  (?)  at  equilibrium  that 


(11) 


In  addition  for  the  case 


Q 

d- 


!!2  = 


dt 


dt 


0  we  get 


(12) 


"3  = 


or  (13) 


n 

n,  =  -  3J  2 


depending  upon  whether  or  not  term  d)  is  neglected.  Under  the  conditions  of 
Interest  Q^>>Q2,  The  ratio  ^^1^2  ®  complicated  function  of  altitude  and 
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Fig.  2.  The  ratio  of  rate  coefficients 
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dep^inds  on  xiy  Hence,  a  direct  calculation  is  impossible.  Wolf  and  Doming  get 
the  vertical  distribution  using  successive  approximation.  Since  the  product 

1; 

Hgn  has  a  scale  height  of  one  half  of  that  for  the  atmosphere,  one  would  expect! 
n^  to  decrease  rapidly  above  the  maximum.  The  decrease  of  above  the  maxi¬ 
mum,  however,  depends  also  on  the  ratio  Q2/Q2  and  hence  on  zenith  angle.  Craig 

-1  '  .  'I  ‘ 

(1950)  gets  a  scale  height  of  approximately  one  half  of  the  atmospheric  scale 

I 

height  for  the  sun  in  the  zenith  and  approximately  the  atmospheric  scale  height 
for  a  zenith  angle  of  6?®,  Equation  (I3)  can  be  reduced  to  the  form 


(14) 

2  q 

”3  ^3  ”3"2'^2  "  *^2  "  ° 

where  (I5) 

b- 

k  =  ^ 

13 

with  (16) 

"^"2  +  V2  \  -  3n2 

n  =  nitrogen  concentration 

”2 

Solving  equation  (14)  for  n^  and  using  the  approximation  that  J2*«12 
there  results 

(17) 

The  equilibrium  time  can  be  estimated  from 
d(n^  +  n-) 

— Tt —  '  %  ■’2  •  ^ V3"i 

Since  n^^nj^  and  the  last  term  on  the  right  of  equation  (18)  is  small,  we  get 
as  a  first  approximation  that 

(19)  n3(t)  =  n3(0)  +  2^  Jgt 

•2 

From  this,  the  time  of  half  restoration  is  estimated  as 

(20)  tf  fir 
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T^la  Is  a 'drastic  sinpliflcation,  however#  since  the  time  oiT  half  restoration 

depends  markedly  on  the  initial  value  of  ny  This  has  been  treat^  thoroughly 

by  Craig  (3.950)1  ami  Taljile  1  tabulates  the  time  constants  estimated  by  him. 

Thie  times  given  (A  t)  apply  to  changes  from  100  n_  to  2.57  n.  ,  2.57  n 

3e  3e 

to  1.50  n  ,  1.50  n  to  1,20  n  ,  1,20  n  to  1,08  n  or  the  reciprocals  of 
3®  3e  3e  3e  3e 

the  foregoing  as  the  case  may  be.  n^^  is  the  photochemical  equilibrium  value. 
Cases  I  and  V  are  for  aenith  angles  of  0°  and  67°  respectively  and  no  dependence 
of  the  oxygen  absorption  coefficient  on  pressure.  Case  II  is  for  aenith  angle 
of  0®  but  viith  a  pressure  dependence  of  p^*^^  for  the  oxygen  absorption  coeffi¬ 
cient. 


TABLE  I 


At  (days) 

Height  (km) 

Case  I 

Case  II 

Case  V 

38 

.15 

.22 

.34 

35 

,60 

.93 

2.1 

32,5 

2.9 

3 

8.9 

30 

13 

11 

52 

27.5 

51 

30 

278 

25 

175 

80 

1,570 

22.5 

650 

184 

4,290  ' 

20 

3,000 

515 

240,000 

17.5 

19,900 

1,470 

15 

84,000 

5,100 

The  case  for  better  agreement  between  the  total  oaone  measurements  and 
the  predictions  of  photochemical  equilibrium  theory  rest  squarly  on  the  fol¬ 
lowing! 

1.  The  dependence  of  the  oxygen  absorption  coefficients  in  the  range 
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1950  A®  to  2150  A®  uppn  pressure.  The  work  of  Hellpern  (19^1 »  1946)  gives  si 
pressure  exponent  of  1.55  at  2144  A®  over  the  range  of  pressure  683  mm  Hg  to 
148  Bin  Hg.  The  pressure  dependence  of  the  absorption  coefficient  In  the  range 
1950  A®  -  2150  A®  dominates  the  production  of  ozone  in  the  range  20  km  to  50  km. 

It  Is  precisely  in  this  range  of  altitude  that  the  majority  of  change  with 

!! 

zenith  angle  of  j'he  sun  occurs  in  the  total  ozone  in  photochemical,  equilibrium. 

2.  An  inclusion  in  the  photochemical  theory  of  more  representative 
temperature  profiles  than  heretofore  have  been  used.  Again  the  region  of  im¬ 
portance  is  20  -  30  km  and  above.  The  necessity  for  this  lies  in  the  extreme 
temperature  dependence  of  the  ratio  of  rate  coefficients  k. 

Some  reasons  for  the  above  statements  will  now  be  given.  With  regard  to 
the  first,  it  la  convenient  to  consider  the  work  of  Craig  (1950)*  Craig  has 
computed  the  total  ozone  expected  at  45°  N  latitude  for  various  solar  zenith 
angles.  The  total  amounts  aret  (Z  *  22®),  ,262  cm;  (Z  =  45®),  .196  cm;  and 
(Z  *  67°),  .185  cm  or  .085  cm,  depending,  respectively,  upon  whether  the  oxygen 
absorption  coefficient  is  or  is  not  pressure  dependent.  Hence,  the  total  ozone 
obtained  is  either  moderately  or  drastically  different  depending  on  whether  the 
absorption  coefficient  is  or  is  not  pressure  dependent. 

In  a  determination  of  the  effect  of  k  on  the  total  ozone,  Craig  used  only 
a  small  seasonal  difference  of  temperature  in  the  stratosphere.  For  example, 
a  temperature  difference  of  only  about  4®C  was  taken  as  representative  for  the 
difference  of  air  temperature  between  the  winter  and  summer  stratosphere  at 
45®  N  latitude.  Furthermore,  an  isothermal  stratosphere  was  assumed  for  both 
winter  and  summer.  This  would  lead  to  a  variance  in  k  of  only  about  20^. 
However,  in  actual  fact,  the  temperature  difference  between  summer  and  winter 
at  45®  N  at  10  mb  is  approximately  20°C.  This  would  produce  a  variation  in  k 
amounting  to  about  320^.  The  two  effects  tend  to  cancel  each  other.  That  is, 
greater  zenith  angles  leading  to  smaller  total  ozone  values  are  accompanied  by 
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larger  values  of  k  leading  to  larger  ozone  concentrations.  Furthermore,  it 
should  be  pointed  out  that  the  observed,  ozone  maximum,  at  least  at  4^^  N,  is 
much  lower  than  that  predicted  by  Craig.  This  lends  more  weight  to  the  seasonal 
variation  of  k  effecting  the  total  ozone  value. 

3.  Measurements  of  the  solar  flux  below  2200  A°.  All  of  the  calcula¬ 
tions  given  by  Craig  are  based  on  an  extrapolated  Aux  in  this  most  important 
part  of  the  solar  spectrum  as  far  as  ozone  production  is  concerned. 

In  spite  of  the  above,  it  appears  that  the  photochemical  theory  alone  will 
not  suffice  to  explain  the  observations. 

In  contrast  to  the  photochemical  production  of  ozone  Just  outlined,  Briner 
(1959)  has  tentatively  suggested  the  following  mechanism  for  the  production  of 
0^  by  mercury  vapor  lamps  by  radiation  in  the  region  2000  A® — >2100  A°. 

(02)2  +  hv— ^0^  +  0 

followed  by  rapid  attachment  of  the  odd  oxygen  atom 
0  +  Og  +  M  +  M 

yielding  art  overall  quantum  efficiency  of  2  ozone  atoms  per  photon  absorbed. 

The  quantity  of  0^^  available  in  the  atmosphere  is  uncertain;  however,  it  has 
been  reported  in  the  telluric  absorption  spectrum  (Goldberg,  1953) •  The  argu¬ 
ment  against  its  existence  in  any  quantity  is  that  its  binding  energy  is  compa¬ 
rable  to  thermal  energies.  The  proportion  of  0^^^  should  be  very  pressure  and 
temperature  dependent.  For  example,  it  constitutes  approximately  50^  of  liquid 
oxygen.  It  is  also  Interesting  along  these  lines  that  the  absolution  coeffi¬ 
cient  for  oxygen  in  the  critical  range  near  2100  A^  is  pressure  dependent  with 
a  varying  exponent,  For  example,  for  the  range  148  -  663  mm  Hg  at  18®C  the 
exponent  is  1.55»  while  from  50  -  130  atmospheres  it  is  2.12. 
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HESULTS 

The  data  presented  and  discussed  in  this  chapter  were  accumulated  solely 
with  the  chemiluminescent  detector.  They  will  be  divided  into  two  parts s 
1.  ground  level  measurements t  2.  vertical  ozone  profiles  obtained  on  balloon 
flights.  Since  the  volume  of  ground  level  data  is  so  large,  only  a  small  por¬ 
tion  of  it  is  presented. 

1.  Ground  level  measurements.  The  surface  measurements  were  made 
in  Minneapolis  at  a  height  of  1  meter  above  the  ground  and  extended  from  Dec¬ 
ember  1,  1961  to  March  1962.  The  most  pronounced  effect  observed  is  the  occur¬ 
ence  of  a  diurnal  variation  on  almost  all  days  but  to  a  varying  degree.  This 
is  shown  in  figures  3  and  4.  The  extent  of  the  diurnal  variation  is  deter¬ 
mined  in  the  main  by  the  magnitude  of  the  surface  level  winds  and  the  sky 
cover.  This  is  understandable  in  terms  of  the  two  types  of  mixing  processes 
which  are  responsible  for  the  arrival  of  ozone  at  the  surface:  mechanical 
mixing,  and  thermal  instability  arising  frim  surface  heating  by  the  sun.  Clear 
nighte  with  low  wind  are  accompanied  by  low  ozone  concentration  near  the  sur¬ 
face.  Clear  days  with  low  wind,  on  the  other  hand,  are  not  necessarily  accom¬ 
panied  by  low  ozone  concentration  because  of  thermal  mixing.  Cloudy  days  and 
nights  are  associated  with  a  small  diurnal  variation  and  an  overall  lower  con¬ 
centration  of  ozone..  The  overall  dependence  of  the  ozone  level  1  meter  above 
the  ground  on  the  surface  winds  for  Minneapolis  throughout  the  month  of  Decem¬ 
ber  is  given  by  the  scatter  diagram  in  figure  5»  Since  the  surface  wind 
measurements  were  those  recorded  by  the  D.S.  Weather  Bureau  at  Wold-ChaniJerlain 
airport  while  the  ozone  measurements  were  made  at  the  University  of  Minnesota, 
a  great  deal  of  scatter  is  to  be  expected.  However,  an  average  relationship 


Ozone  measurements  1  m  above  the  grotmd  at  Minneapolis  on  Dec.  23*  24  and  25*  1961 
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between  surface  wind  and  ozone  concentration  is  clear.  None  of  the  ideas  in 
this  section  make  sense  without  the  existence  of  an  ozone  sink  at  the  earth's 
surface,  which  is  continually  removing  ozone  from  the  atmosphere.  The  large 
diurnal  variation  observe  near  the  surface  prompted  measuremerits  of  the  verti- 

i  ' 

cal  distribution  by  mearis  of  a  tethered  balloon  which  could  be  raised  and 
lowered  through  the  lower  lasrer.  This  work  was  done  and  reported  by  Kroenlng 
and  Ney  (1962).  For  convenience,  it  is  reproduced  in  figure  6.  The  graph 
shows  three  vertical  distributions  of  ozone  spaced  in  time  shortly  after  sun¬ 
set  which  were  obtained  at  the  Anoka  County  airport  on  June  20,  196I.  It  was  a 
clear  night  with  low  surface  iflnds.  The  vertical  temperature  distribution  was 
monitored  at  the  same  time.  The  model  .adopted  from  this  data  is  that  idien  con¬ 
ditions  allow  a  temperature  inversion  to  form  near  the  surface,  the  process  of 
decomposition  removes  the  ozone  from  the  layer  near  to  the  surface.  This 
provides  a  simple  and  direct  measurement  of  the  sink  strength.  From  the  data 
in  figure  6,  Kroening  and  Ney  (1962)  have  deduced  a  sink  strength  of  6  x  10^® 
atoms /om^  sec  by  merely  noting  the  quantity  of  ozone  removed  and  the  time 
required  for  this  to  occur.  Although  this  value  of  the  sink  strength  is  really 
an  average  value  over  a  specified  time  interval,  it  is  difficult  to  allow  a 
variance  of  more  than  a  factor  of  2  in  the  above  estimate. 

After  establishing  the  value  of  the  sink  strength  as  in  the  above,  it  is 
possible  to  determine  the  kinetic-theory  properties  associated  with  this  sink. 
Since  the  eddy  diffusivity  in  the  atmosphere  near  the  ground  is  much  larger 
than  the  kinetic  theory  value,  the  vertical  ozone  distribution  is  determined 
completely  by  eddy. diffusion.  Tet,  in  the  last  few  mean  free  paths  kinetic- 
theory  phenomena  must  take  over.  Thus,  the  value  of  the  sink  strength  can  be 
used  to  determine  the  number  of  collisions  required  before  decomposition  occurs. 
The  appropriate  equation  Isi 
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Fig.  6.  The  ozone  concentration  and  air  temperature  vs  altitude 
in  the  lower  surface  layers  shortly  after  sunset  on 
June  20,  1961. 
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(21)  l/4iuvf  =  6  X  10^°  S.t|ms 
^  ■  cnrsec 

■where  is  the  number  density  of  ozone  molecules,  v  is  the  mean  molecular 
velocity  and  f  is  the  fraction  of  collisions  with  the  surface  leading  to  decom¬ 
position  of  an  ozone  molecule •  Putting  in  the  appropriate  numbers 
n^  =  6  X  10^^/cm^ 

V  *  3.6  X  10**  cm/seo 

f  is  determined  to  be  approximately  lO"-^,  Hence,  only  one  collision  in  100,000 
is  effective  in  decomposing  an  ozone  atom  and  removing  it  from  the  atmosphere. 
This  is  characteristic  of  many  chemical  reactions  in  which  aotivatii^jn  energies 
can  be  overcome  only  by  molecules  on  the  tall  of  the  maxwelllan  destrlbution  in 
order  for  the  reaction  to  proceed.  The  value  of  f  =  10"^  falls  within  the 
range  of  values  measured  for  ordinary  laboratory  materials.  For  example,  f 
for  black  rubber  is  2  x  10"**,  while  for  mayon  tubing  f  »  2  x  10"®.  Mayon 
tubing  is  so  resistant  to  the  decomposition  of  ozone  that  it  can  be  used  in 
ozone  plumbing  with  negligible  loss. 

With  the  value  of  f  at  hand,  it  is  now  possible  to  back  track  and  make  the 
following  point.  In  the  above  discussion,  it  was  clearly  demonstrated  that  the 
decomposition  rate  of  ozone  near  the  earth's  surface  is  diffusion  controlled. 
The  mathematical  formulation  of  this  is  as  follows: 

(aa) 

Where  P  is  the  vertical  diffusive  flux  of  ozone,  D  is  the  vertical  eddy  dlf- 

dn- 

fusion  coefficient  and  is  the  vertical  gradient  of  the  ozone  concentration. 
If  the  above  vertical  flux  is  equated  to  the  strength  of  the  sink,  we  obtain 
dn  (z) 

- 1/**  "3(0>^ 

This  equation  relates  the  vertical  eddy  diffuslvity,  the  vertical  ozone  con- 
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centration  gradient  in  a  simple  manner  once  f  is  knovm.  It  can  Tae  applied  where 
vertical  diffusion  Is ihe  sole  transfer  proooess  and  only  in  that  part  of '  the  atmos¬ 
phere  where  ozone  is  a  conservative  constituent.  Under  steady  conditions,  we 
would  then  have  the  production  at  high  levels  balanced  by  the  destruction  at 
the  surface.  In  this  state,  a  quantity  of  ozone  would  be  stored  in  the  atmos¬ 
phere  which  depended  solely  on  the  production  at  high  levels,  the  eddy  dif- 
fusivity  and  the  product  vf. 

Equation  (23),  although  written  for  the  steady  state,  lends  some  insight 
into  the  non-steady  state  process  near  the  surface.  If  conditions  at  some 
locality  are  such  that  ^  does  not  change  appreciably  over  a  certain  time  inter¬ 
val  but  the  sink  strength  changes  in  accordance  with  n^(0),  it  is  possible  to 
generate  high  ozone  concentrations  in  the  troposphere  near  the  ground.  If  there 
la  a  change  of  the  steady  state  conditions  to  which  (23)  applies  to  the  case 
where  the  sink  effectively  vanishes,  as  would  be  in  the  abrupt  formation  of  an 
inversion  near  the  ground  or  stagnant  air  masses,  the  ozone  concentration 
Immediately  begins  to  increase  above  the  inversion  or  stagnant  air  mass.  If 
these  conditions  persisted  long  enough,  the  ozone  would  back  up,  so  to  speak, 
all  the  way  through  the  atmosphere.  However,  the  majority  of  these  situations 
are  short  termed,  and,  furthermore,  horizontal  transport  will  wash  out  the 
effect  by  horizontal  diffusion  and  adveotlon.  Yet,  there  should  be  an  average 
effect  of  this  kind  for  the  whole  atmosphere,  and,  certainly,  it  should  be 
visible  in  the  individual  oases  as  the  following  calculation  shows: 

(24)  An-(z  J  » 

3  o 

Substituting  the  appropriate  number,  equation  (24)  yields 

6  X  10^*^  Bt|jns  ^  lo^sec 

- - 2E-S22 - -  =  6  x  10^°  ^ 

lO^cm  ®® 

This  calculation  is  consistent  with  the  downward  flux  and  the  vertical  eddy 
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iji  dif fusivlty  in  the  troposphere.  The  build  up  is  considere^^  to  take  place  in 

”\the  time  interval  of  a  day  and  within  a  layer  1  km  thick.  J  it  was  also  assumed 

1  -  1 

that  the  ozone  mixing  ratio  increases  by  10^  per  km.  The  data  shown  in  figures 

'•  ‘-t  '!■  I,  '■ 

3  and  4  are  consistent  with  such  an  explanation.  The  graphs  irioiude  days 
(12/20/62  and  12/24/62)  irtien  the  resultant  wind  at  the  sutfaea  is  very  small. 

This  is  preceded  and  followed  by  days  in  which  there  was  good  movement.  The 
days  after  the  day,  of  low  resultant  wind  are  accompanied  by  higher  than  normal 
ozone  at  ground  level.  This  trend  is  the  same  for  9  out  of  I3  instances  on 
record.  Hence,  it  is  clearly  brought  out  once  again  that  the  ozone  content  of 
the  atmosphere  is  dependent  on  surface  conditions  via  a  diffusion  controlled 
process  in  the  layers  near  to  the  earth's  surface.  Still  there  is  another 
point  to  be  made  regarding  the  sink  strength,  and  it  is  probably  the  most  impor¬ 
tant.  One  is  forced  lo  consider  the  variation  of  ^  in  the  sink  strength. 

1/2 

Although  V  is  proportional  to  T  '  where  T  is  the  aniblent  surface  temperature 
in  degrees  Kelvin,  the  factor  f  will  have  an  exponential  temperature  dependence 
and  will  certainly  depend  on  the  nature  of  the  surface.  But  it  is  not  obvious 
how  one  can  handle  this  last  variance.  For  the  discussion,  it  will  be  assumed 
that  there  is  some  substance,  perhaps  water,  which  may  be  considered  as  repre-  1 

sentativs  for  the  whole  surface.  That  this  assumption  is  valid  might  be  indi¬ 
cated  by  a  northern-southern  hemisphere  total  ozone  difference  since  the 
northern  hemisphere  contains  much  more  land  mass.  It  will  also  be  assumed  that 
the  sink  strength  determined  at  Minneapolis  from  figure  6  is  representative  of 

■j 

the  activation  energy  for  decomposition  at  that  temperature  and  that  this  i 

activation  energy  is  the  same  for  the  whole  northern  hemisphere.  Thus,  ! 

(25)  =  lO"'^  I 

if  there  is  no  steric  hindrance  in  the  decomposition  process.  is  the  energy 
of  activation,  T  is  the  absolute  temperature  and  k  is  the  Boltzmann  factor. 

With  T  =  300®K  the  solution  of  (25)  is  E^  =  .3  ev.  With  E^  in  hand,  we  can  con- 
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struct  thei  temperature  dependence  of  f.  This  is  given  In  figure  7.  Hence, 

1)  ",  .  , 
finally!;  have  a  mean?  of  treating  the  variability  of  the  sink  strength. 

Assuming  an  activation  energy  of  .3  ev,  equation  (25)  predicts  a  variation  of  a 

factor  of  300  for  a  temperature  change  of  200®K  to  300®K,  Clearly,  then,  the 

sink  strength  must  be  Included  In  any  equilibrium  theory  of  atmospheric  ossone. 

We  will  now  proceed  to  estimate  the  effect  of  the  sink  on  the  total  ozone. 

Rewriting  equation  (23)  so  as  to  allow  for  variation  In  air  density  results  in 


(26)  da  "  n2(o,X)vf 

where  D  =  -2 

n 

If  1)  source  and  sink  are  In  equilibrium 

2)  horizontal  transfer  Is  unimportant 

3)  furthermore  molecular  diffusion  Is  unimportant 

4)  we  consider  only  that  part  of  the  atmosphere  in  which  ozone  is  a 
conservative  constituent,  then  Integration  of  (26)  yields 

(27)  nJz,\)  =  n3(0,X)e“®''^  +  e"®/^  \ 

3  3  'oKD 

for  the  vertical  ozone  distribution.  Integrating  (27)  between  z  =  o  and  z  =  h 

2 

gives  the  total  ozone  In  a  1  cm  column  of  height  h  as 


h/H 

(28)  N-(h,\)  *  n-(o,X.)H(l-e"  )  + 

J  J  '0 


h 


(  n  (o,X)^ 


dz 


In  the  above  n  » 
X  s 
H  = 
D  = 
n2(z,X)  = 
N^(h,X)  = 


•<21 /H 

n^e”  '  hydrostatic  equation 
latitude 

atmospheric  scale  height 

vertical  eddy  dlfl!Usivlty 

ozone  number  density 

total  ozone  between  z  =  o  and  z  ^  h 


(i<FA»ann«i|)(r«i 


-  ..^\ 


k 

1 


Taking  h  =  kU 


ii9)  .  \  N^(ilfl,X)  l;n^(o,X)H+ j 


n  (o,X)vf 


ds  dz 


If  it  is  now  assumed  that  the  source  strength  and  D  are  independent  of  latitude* 
then  the  last  thrm  on  the  rlghi  of  (29)  is  a  constant.  The  first  term  then 


represents  the 


ozone  stored  in  the  atmosphere  to  compensate  for  the  variation 


:  in  f .  This  te^,)in  will  be  dependent  on  latitude. 


n„(o,X) 


(30)  =  N^C^H.X)  -  N3(4H.0)  =  H[nj(o,X)-nj(o,0)]=  Hn^Co.O)  -  ij 

The  latitudinal  dependence  of  this  tern  is  shown  in  figure  8*  curve  B*  eemputed 
from  London's  (1957)  mean  temperatures  and  under  the  assumption  that  the  pro> 
duct  n3(o,X)  f  is  a  constant.  If  curve  B  is  weighted  according  to  the  yearly 
mean  insolation  as  given  by  Fritz  (1951) •  curve  C  is  the  result.  In  a  similar 
manner  one  can  obtain  the  seasonal  variation 


The  ratio of  the  seasonal  variatim  to  the  mean  latitudinal  variation 
is  given  in  figure  8,  curve  A.  Suffice  it  to  say,  however,  that  the  magnitude 
of  these  terms  is  small  compared  to  the  total  ozone  observed.  This  is  shown  in 
figure  9  which  gives  the  observed  yearly  mean  latitudinal  variation  while  the 
solid  curve  gives  the  expected  variation  from  the  above  analysis,  assuming  the 
total  ozone  at  the  equator  is  240  mllli-atmo-cm  and  that  n3(o,o)  =  4x10  /cim'^. 
The  seasonal  variation,  while  small,  is  in  phase  with  the  temperature.  Althou§^ 
the  above  considerations  appear  unimportant  regarding  total  ozone,  they  still 
predict  a  latitudinal  variation  in  the  ozone  concentration  in  the  troposphere 
of  a  factor  of  about  six  from  the  equator  to  the  pole. 

Hence,  we  must  look  for  variations  occuring  in  the  second  term  on  the  right 
side  of  equation  (29)  ♦  This  term  can  be  written  under  the  previous  assumption 


1 


as: 


(32) 


n^(otX-)yf 


-39- 
Hs 


4H  Pi 
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-s/H 


ds 


dz 


This  i^nia  may  ^  reduced  to.  the  approximate  form 

(33)  _2__ -  X  2i_ 

«  .  4,  D 

where  D  Is  the  average,  vertical  diffusion  coefficient  over  the  range  z  =  o  to 

V' 

z  a  4H,  Cutting  off  the  integral  at  z  =  4H  places  the  source  at  this  altitude 
(about  2?  km)  and  neglects  the  ozone  above  this  altitude.  Further,  the  observed 
yearly  mean  total  ozone  values  would  be  accounted  for  by  this  term  in  a  craa- 
pletely  diffusion  controlled  theory.  The  yearly  mean  value  of  5  required  at 
45®  N  is  given,  for  example,  by 


m,(o,X)3H^ 

(34)  -  - -  =  2.68x10^® 

4  D 


X  3.5 


Substituting  the  previously  used  sink  strength  and  H  *  7.5  x  10^  cm 

2.68  X  10“  X  3.5 

The  mean  values  of  D  required  as  a  function  of  latitude  are  given  in  Table  II. 
The  magnitude  of  D  depends  on  the  reciprocal  of  the  yearly  mean  total  ozone. 
Furthermore,  the  seasonal  total  ozone  might  similarly  be  explained  by  the 
seasonal  variation  in  5  given  in  Table  II.  According  to  this  equillb?d.um 
theory,  <equatlon  (26)  would  also  yield  the  vertical  eddy  dlffuslvity  ffom  a 
measurement  of  the  vertical  profile.  The  eddy  diffusivlty  would  then  be  pro¬ 
portional  to  the  reciprocal  of  the  gradient  of  the  mixing  ratio,  l.e. 


(35)  D  l/4n(o.A)vf 


dz 


dz 


l/4n(o.X)^ 


3rT^(t/L  Mg  dz 


If  n^  is  expressed  in  atoras/cm^  and  AT  In  ®K,  the  above  becomes  for  5  tan  inter¬ 
vals 
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(sec 


6  X  10^° _ V  ■' 

An3x2ri.0-^+n2x(^1.23xl0’^^^l  +  0J 

J  -  ,  _  .,  .  ■ 

TABLE  n  ~ 


Coniputed  Mean  Vertical  Difflision  Coefficient  end  Seasdnal  Variation  Thereof 


Lati¬ 

tude 

0 

10 

20 

30 

N  (ffill- 
11-atmo- 

150 

170, 

200 

275 

cm) 


40 

50 

6,0 

70 

80 

325 

350 

350 

350 

320 

6  (X)  cm^  2.5x10^  2.2x10^  1.9x10^  1.4x10^  1.2x10^  1.1x10^  1.1x10^  1.1x10^1.2300^ 
per  sec 


AD  q) 

5  (X) 


0  0 


0  .1  .18  .28  .28  .30  .4 


The  values  computed  for  flight  600  In  figure  11  are  tabulated  In  Table  m. 
5  km  smoothed  Inteinrals  were  used  In  the  calculations. 


_ _  TABLE  in _ 

_ Computed  Vertical  Eddy  Dlffuslvlty  for  Flight  600 

_ 0"3km _ 5-lOkm  10-15km  15-20km  20-25km 

^5km  {mcI  x.lO^  6.3  X  10^  1.2  x  10^ 


The  values  obtained  for  D  and  are  not  unreasonable  in  view  of  the 
estimates  of  this  quantity  by  Lettau  (1951).  Hence,  the  above  calculations 
demonstrate  the  possibility  that  the  majority  of  the  ozone  distribution  may  be 
diffusion  controlled  and  that  horizontal  diffusion  and  advectlon  may  be  unimpor¬ 
tant  for  the  ’gross  features . 

By  far,  the  most  significant  feature  of  the  ozone  measurements  Is  their 
coirespondenee  with  the  temperature  profile.  This  starts  at  the  ground  where  it 
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is  noticed  that  a  decrease  in  ozone  acoompahids  a  temperature  inversion  near 

yj  ' 

the  surface.  Moreover*  the  observed  mixing  ratio  in  the  tropospheie’  depends 
on  the  temperaiure  lapse  rate.  The  most  pronounced  variation  in  the  ozone 
mixing  ratio  occurs  at  the  tropopause.  A  well  defined  tropopause  (abiropt  r 
change  in  the  temperature  lapse  rate)  is  accompanied  by  a  very  sha^  increase 
in  the  ozone  mixing  ratio.  Likewise*  where  the  temperature  lapse  rate  changes 
gradually*  the  ozone  indicates  a  poorly  defined  boundary  between  troposphere 
and  stratosphere.  Hence,  the  diffusive  transport  must  be  controlled  to  a  large 
extent  by  the  temperature  distribution,  and  it  is  the  semi-impermeable  surface 
termed  the  tropopause  which  hinders  any  large  quantity  of  ozone  from  entering 
the  troposphere.  This  diffusion  controlled  process  is  advocated  up  to  the  level 
where  photochemical  time  constants  become  comparable  to  the  residence  times 
consistent  with  the  diffusive  process.  This  level  is  probably  in  the  neigh- 
bdrhood  of  25  km  in  the  vicinity  of  the  observed  ozone  maximum  which  is  lower 
than  the  photochemical  maximum.  This  Interplay  between  diffusive  and  photo¬ 
chemical  processes  is  suggested  as  the  cause  for  the  observed  ozone  maximum 
being  lower  than  the  photocho&ical  maximum.  Furthermore*  the  concept  of  a  net 
ozone  production  rate  independent  of  the  solar  zenith  angle  and  latitude  is  not 
inconsistent  with  the  above  ideas.  The  calculations  of  Craig  (1950)  showi 

1.  a  maximum  ozone  concentration  at  30  km  for  small  zenith  angles; 

2.  a  rapid  decrease  of  the  ozone  mixing  ratio  above  the  maximum  with 
a  scale  height  of  roughly  half  of  the  atmospheric  scale  height.  The  second  of 
these  is  characteristic  of  any  photochemical  equilibrium  theory.  The  first  is 
dependent  upon  assumptions  regarding  the  oxygen  absorption  process.  However* 
neither  of  the  above  are  verified  by  the  observations  in  this  work.  With  a 
large  diffusion  coefficient  and  with  the  rapid  decrease  of  the  ozone  mixing 
ratio  abcve  the  maximum  required  by  photochemical  equilibrium*  ozone  would  be 
transported  upward  where  it  would  be  photochemicaUy  destroyed.  This  process 
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would  lowei!r  the  maxinun  and  increase  the  oaotle  scale  height- above  the  maximum. 
For  large  solar  zenith  angles,  however,  the  photochemical  scale  height  iils 
approxtoately  that  of  the  atmosphere,  greatly  reducing  or  preventing  upward 
transport.  Hence,  it  is  possible  that  large  solar  zenith  angles  yield  no  less 
ozone  above  the  maximum  than  small  zenith  angles  unier  identical  tenpiprature 
and  diffusion  conditions.  ''  To  reverse  the  effect  under  identical  diffusion  con^ 
ditions  would  be  possible  since,  as  was  outlined  earlier  in  the  section  on 
photochemical  equilibrium,  the  k  factor  is  very  temperature  dependent  and  is 
such  as  to  increase  the  ozone  concentration  at  high  levels  ih  the  winter  or 
early  spring  and  to  deorease  it  in  the  late  summer.  The  observations  shown  in 
figures  10,  11,  12  and  13  support  this.  Moreover,  if  the  absorption  coefficient 
of  oxygen  is  pressure  dependent,  the  above  argument  would  lead  to  a  broad 
maximum  for  large  solar  zenith  angles  and  a  more  pronounced  maximum  for  small 
zenith  angles.  This  also  is  suggested  by  the  data  in  figures  10,  11,  12,  and 
14.  It  is  reasonable  to  expect,  then,  that  the  observed  maximum  exists  at  a 
level  where  the  photochemioal  time  constant  and  the  residence  time  for  the 
transport  process  are  the  same.  This  level  may  change  with  the  season.  The 
data  in  figures  10,  11,  12,  and  13  suggest  that  this  level  is  higher  in  the 
spring  (broad  maximum)  than  it  is  in  the  summer  (sharp  maximum).  This  is  cion> 
sistent  with  a  larger  photochemical  time  constant  at  large  solar  zenith  angles. 
However,  it  is  mandatory  that  more  observations  be  made  to  test  these  ideas. 

The  vertical  profiles  of  the  atmospheric  ozone  distribution  shown  in 
figures  10,  11,  12,  and  13  are  unique  in  two  respects! 

1.  They  are  among  the  first  of  such  measurements  which  yield  struo- 
ture  in  the  vertical  ozone  distribution  heretofore  unresolved  while  at  the  same 
time  providing  a  measure  of  the  integrated  ozone  between  ground  level  and  bal¬ 
loon  altitudes. 

2.  They  furnish  measurements  of  ozone  in  a  manner  which  is  com- 
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pletely  different  itom  other  techniques. 

The  chemiluninescerit  ozone  detector  will  easily  resolve  structure  in  the 
vertical  distribution  viiieh  cannot  be  shown  in  figures  10,  11^  12,  and  13. 

This  teatu|jre  is  a  prerequlslte^th  the  study  of  advective  processes.  It  has  beeri' 
shown  (Kroening  and  Ney,  1962)  that  ozone  structures  with  vertical  dimensions 
of  the  order  of  1  km  exist  in  the  troposphere.  These  thin  stable  laminae, 
termed  ozone  'rivers,'  are  attributed  to  the  vertical  advection  of  air  between 
the  stratosphere  and  the  troposphere  and  thus  retain  an  ozone  concentration 
characteristic  of  their  origin.  These  ozone  rivers  are  shown  in  figure  10  and 

11  and  appear  to  be  characterlstle  of  the  late  summer  soundings  (July  and  « 

1  ! 

August).  The  distinguishing  feature  between  an  ozone  river  and  a  diffusive 
maximum  has  been  assumed  to  be  in  the  nature  of  the  mixing  ratio  above  the 
maximum.  In  the  formation  of  a  diffusive  maximum,  the  mixing  ratio  cannot 
decrease  above  the  maximum.  In  the  ease  of  advection,  this  criterion  need  not 
be  satisfied  and  most  likely  will  not,  For  example,  the  maxima  at  1^5  mb  in 
figure  13  and  I65  mb  in  figure  10  are  indicative  of  an  advective  process, 
that  is,  an  ozone  river.  In  neither  of  these  cases  is  it  clear  what  fraction 
of  the  advection  is  vertical  and  what  fraction  is  horizontal.  Most  likely  it  is 
a  combination  of  both.  The  temperature  profile  associated  with  these  rivers, 
however,  suggests  downward  transport  as  the  major  advection.  Certainly,  soun¬ 
dings  on  a  synoptic  basis  or  a  vehicle  navigating  on  these  rivers  would  provide 
the  answer.  It  has  been  shown  (Ney  and  Kroening,  1961)  that  these  ozone  rivers 
can  only  persist  in.  the  case  of  continued  descending  motion  and  it  was  this  idea 

i; 

that  led  to  their  naming.  Any  ascending  lamina  will  be  characterized  by  a  low 
ozone  concentration  and  much  larger  vertical  dimensions. 

The  three  balloon  flights  (593i  600,  596),  besides  providing  the  vertical 
ozone  distribution,  were  maintained  at  high  altitude  through  sunrise  in  order 
to  determine  the  diurnal  ozone  variation  at  high  altitude.  The  resutls  obtained 


13.  Vertical  ozone  dlstribatlon  and  taaperature  profile  over  Minneapolis  on  April  28,  1962. 
Carve  A  is  the  partial  pressure  of  0_  in  pnb.  Curve  B  is  the  ozone  density  in  pg/nP 
with  sane  scale  as  A. 


ruled  out.  There  is  a  tendency  for  a  dovmtmH  flow  of  air  around  the  halloon 
at  night  and  the  opposite  eiroulation  after  sunrise.  This  could  explain  the 
higher  readings  after  sunrise.  Certainly,  the  diurnal  variation  is  less  than 
3056,  More  refined  experiments  are  evidently  needed, 

-  i 
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PRCDDCTION  OF  OZONE 

The  source  of  ozone  as  outlined  earlier  in  this  vnriting  was  taken  as  that 
oocsuring  in  the  high  atuosphere  from  solar  energy.  It  appears,  however >  that 
this  cannot  be  taken  as  the  only  significant  source  (Kroening  and  Ney). 
Electrical  discharge  phenomena  apparently  involve  enough  energy  dissipation  with 
a  high  enough  efficiency  to  compensate  for  the  ozone  removed  from  the  atmosphere 
at  the  ground.  Cauer  (1951)  gives  strong  evidence  for  the  production  of  ozone 
in  regions  of  high  electric  fields  without  the  spark  discharge.  Regener  (19^) 
has  reported  observations  of  high  ozone  levels  under  thunderstorm  oonditidns  at 
cloud  levels  in  the  troposphere.  Kroening  and  Ney  (1962)  have  measured  an 
efficiency  factor  for  ozone  production  by  a  radioactive  source  and  arrive 
at  the  value  .03  atoms  of  ozone  per  electron  volt  dissipated  in  the  air.  This 
is  about  one  ozone  atom  for  every  ion  pair.  A  nearly  identical  efficiency 
factor  Is  obtained  in  the  laboratory  controlled  eleotrieal  discharge  in  air. 

It  is  thought  that  this  is  more  than  coincidental  and  suggests  the  conclusion 
that  this  efficiency  factor  applies  to  any  very  energetic  process.  In  light  of 
the  possibility  of  the  production  of  large  amounts  of  ozone  by  electrical  dis¬ 
charges.  experiments  should  undoubtedly  be  carried  out  to  check  this.  Sreasure- 
ments  around  thunderstorms  should  be  made  aloft.  Surface  measurements  arf»  like¬ 
ly  to  be  obscured.  However,  it  should  be  mentioned  that  during  a  large  thunder¬ 
storm  occuring  on  5/22/62  and  5/23/62,  a  very  large  increase  in  ozone  arriving- 
at  the  surface  was  recorded  at  Minneapolis.  The  ozone  reached  a  level  sllghia.y 
in  access  of  10  pphm  by  volume  while  the  average  for  an  18-hour  periocf  was 
6  pphm,  almost  twice  the  mean  surface  level.  In  view  of  the  data  available  at 
present  concerning  the  production  of  ozone  in  electric  dlsohai'ges,  it  appears 
that  it  cannot  be  neglected. 


SDMMABI 


1.  A  ohemilunlnesoent  ozone  detector  is  described  which  allows 

i' 

detailed  ineasarements  of  the  atmospheric  ozone  concentration. 

2.  There  is  a  diurnal  variation  of  the  atmospheric  ozone  concentra¬ 
tion  1  m  above  the  ground  on  almost  all  days  but  to  a  varsring  degree. 

3.  An  average  relationship  between  the  surface  level  wind  speed  and 
the  ozone  concentration  1  m  above  the  ground  at  Minneapolis  is  given. 

4.  The  strength  of  the  ozone  sink  at  Minneapolis  is  given  and  from 
this  the  general  properties  of  the  ground  level  sink  are  deduced.  A  definite 
latitudinal  variation  of  the  tropospheric  ozone  concentration  results  with  more 
ozone  at  higher  latitudes.  The  effect  of  the  ground  level  sink  on  the  total 
ozone  value,  however,  is  shown  to  be  negligible  if  an  allowance  is  madb  for  the 
solar  insolation  available  at  the  top  of  the  atmosphere. 

5<  Theoretical  and  experimental  evidence  for  the  formation  of  a 
low  level  ozone  maximum  in  the  troposphere  is  given. 

6,  A  completeljf'  diffusion-controlled  model  of  the  vertical  ozone  dis¬ 
tribution  is  presented  and  the  appropriate  values  for  the  vertical  eddy  diffttsi- 
vity  are  given, 

7,  The  existence  of  thin  stable  lamina  of  ozone,  termed  ozone 
'rivers,'  la  demonstrated.  Those  rivers  descend  through  the  atmosphere  and  re¬ 
tain  an  ozone  concentration  characteristic  of  their  origin. 

8,  It  is  pointed  out  that  the  production  of  ozone  by  electrical  dis¬ 
charges  within  the  atmosphere  cannot  be  ruled  as  insignificant. 
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Procftdupe  for  the  Chemloal  Detennlnation  of  ,  i 

.  ,■  ,  'i 

Ozone  Conoentrationa 

Reagents! 

Buffer  solutions  -  1.8  g  of  NSg  HPOj^  y  +  1,7  g  of  KHgPOK  In  1  liter  of 
water . 

.OOOINOO  N  or  .OOlOON  KIG  -  .0214  g  or  .214  g  or  KIO.  in  1  liter  of 
water.  3 

2  N  sulforic  acid. 

KI  free  from  iodate. 

AmperoBietrio  probe,  magnetic  stirrer,  5  ml  buret  graduated,  Flovmieter, 
bubbler  bottles. 

Procedure I 

Place  70  ml  of  buffer  solution,  1  g  of  KI  and  ^  ml  of  NsgSgO^  (measured 
with  5  ml  pipette)  solution  in  each  of  the  two  bubbler  bottles.  Pass  ozone 
through  at  the  rate  of  500  co  per  minute  for  about  20  minutes.  Add  10  ml 
ZN-HgSO^i^.  Titrate  with  standardized  KIO^  solution  from  buret  using  an  ampero- 
metrlc  probe  made  of  silver  and  gold  (anode  and  cathode  respectively). 

Caution  should  be  taken  to  cleanse  all  apparatus  before  the  sample  Is  taken. 
Distilled  water  should  be  used  for  all  solutions,  A  dry  run  is  generally  made 
first  (without  ozone)  to  test  the  stability  of JNSgSgO^  during  the  sampling  time. 
Care  must  be  taken  not  to  use  all  of  the  Na^SgO^  during  the  ozone  determination 
(jr  a  rerun  with  more  NsgSgO^  or  less  time  is  required,  A  blank  a  day  is  a  good 
idea  to  keep  the  decomposition  of  the  solution  in  view.  The  bubblers 

used  should  be  of  the  large  hole  type.  Fine  hole  bubbler  bottles  decompose 


ozone.  The  plumbing  should  be  ell  glass  and  plastio.  No  stopcock  greasy  should 


be  used  on  fittings.  In  aU  of  the  calibrations*  tank  oxygen  was  used. 
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Surface  Temperature 

—  London’s  data 

Latitude 

Yearly  Mean  Temp 
(°G) 

Winter  Temp 

(®c) 

Summer  Temo 

(°c) 

0-10 

26.8 

27.3 

27.3 

0-20 

26.8 

24.7 

27.9 

20  -  30 

23.7 

19.4 

26.9 

o 

i 

o 

17.4 

10.5 

23.2 

40-50 

9.5 

.4 

18,3 

50  -  6o 

2.6 

-8.5 

13.3 

6o  -  70 

-5.3 

-18,9  - 

9.6 

70  -  80 

-11.3 

-25,1 

3.4 

80  -  90 

-16.4 

-30.2 

.4 
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i  Introduction 

The  equillbrluni  theory  of  small  atmospheric  ions  and  the  technique  of 
measuring  them  has  been  outlined  elsewhere  (Kroening,  1959) <  Ihe  theory  in¬ 
volves  the  following  equation! 

(1)  ^  =  q  -  o<n2  -  I  p(R)  N  (R)  n 

I 

where  n  is' the  number  density  of  small  ionsi  N  is  the  number  density  of  large 
nuclei,  charged  or  uncharged,  q  is  the  production  rate  of  small  ions,  is  the 
Thomson  volume  recombination  coefficient  and  p  is  the  related!  coefficient  for 
recombination  between  small  ions  and  nuclei.  The  production  rate  is  usually 
attributed  to  cosmic  rays.  The  last  term  on  the  right  side  of  the  equation  (l) 
has  been  neglected  or  deemed  unimportant  in  the  equilibrium  theory  above  the 
esEOhange  layer.  It  has  been  neglected  largely  because  an  understandings  of  the 
SDiall  ion  ^  large  nuolei  interaction  is  laoking  except  at  ground  levels,  and 
even  here  the  value  of  P  is  not  a  well  defined  quantity.  Eor  example,  the 
values  of  p  reported  vary  from  .il’  x  lO'^om^/sec  up  to  8  x  lO'^cn^/seo.  An 
average  value  of  2  x  10''”om'’/sec  at  H  T  P  is  normally  adopted,  Purthennore,  it 
is  the  opinion  of  the  author  that  any  conductivity  measurements  or  ion  density 
measurements  made  at  heights  greater  than  25  km  in  the  daytime  are  subject  to 
serious  error  unless  special  precautions  have  been  taken  to  eliminate  photo- 
electron  ejection  from  the  collector.  This  effect  has  been  observed  on  every 
daytime  flight  made  by  the  author  and  if  not  eliminated  obscures  the  effect  of 
large  nuclei  (Kroening,  i960).  Mention  of  this  in  the  literature  has  been  made 
only  by  the  author.  Yet,  daytime  measurements  which  are  likely  in  error  are 


n 


taken  as  representative  of  ambient,  conditions*  Moreover*  tbe  oonducitivity  „ 
measurements  are  not  djjeotly  applicable  to  the  equili|brlum  'ihepi7  Ance  these 
measurements  have  the  ipnie  mobility  folded  into  them.  Thb  ion  density  measure¬ 
ments  reported  in  this  writing  and  elsewhere  were  made  at  night  to  eliminate 
these  i[iiroUems.  Consequently*  it  is  advocated  that  they  are  representative  of 
ambient  conditions  in  the  high  atmosphere.  Daytime  measurements  of  the  conduc¬ 
tivity  at  high  altitude  which  are  explained  solely  by  small  ion  recombination 

li 

must  be  Judged  as  fictitious*  since  all  of  the  authors  measurements*  except 
one*  show  a  deficit  in  the  small  ion  density  at  high  altitudes.  It  iij  felt 

that  this  can  only  be  explained  by  the  presence  of  Aitken  nuclei  at  high  alti- 

1 

tude. 
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THEOET  AI1D,,MEASUEE»BNTS 


The  snail  ion  reoonbination  process  has  heeh  studied  oonsiderably  and  the 

V, 

recombination  eoeffioient  governing  it  has  been  given  bjr  Thonaon*  Very  little 
is  known  coneeming  the  small  ion  -  large  nuolei  recombination  proeess>  Eenoe* 
in  this  part  of  the  disserta^on*  an  attempt  wili  be  made  to  describe  the  ijspor« 
tanoe  of  the  last  tern  on  the  right  of  equation  (1) •  eKpeoially  at  high  alti> 

tudes. 

Measurenents  of  q  for  the  appropriate  latitude  (45®  N)  as  well  as  a  graph 
of  a  have  been  given  elsewhere  (Kroening,  1959) .  In  an  attempt  to  estimate  ^ 
a  simple  nethod  of  arriving  at  the  value  of  a  will  be  presented.  The  attrae- 
tive  force  between  positive  and  negative  singly  charged  ions  is  e  /r  >  whereas 
the  potential  between  the  two  is  e^/r.  For  the  two  ions  to  recombine  requires 
that  the  relative  kinetic  energy  is  such  as  to  'allow  a  closed  orbit  of  wie  about 


the  other.  In  a  gas  at  temperature  T,  this  requires 

e^  ^  IkT 
r  2 


(2) 


Hence f  they  must  approach  at  least  as  close  as 


(3) 


r  ® 

®  3kT 


Furtheimorey  during  the  tine  that  the  two  ions  satisfy  equation  (2),  a  collision 
of  either  of  the  two  ions  with  an  air  molecule  which  reduces  the  relative  kine¬ 
tic  energy  by  an  amount  kT  must  occur.  Then  recombination  is  effected.  These 
are  the  fundamental  assumptions  made  by  Thomson.  If  the  probability  of  the 

collision  of  the  ions  with  the  neutral  air  molecules  while  r  ^  r^:  are  s,,.  and  s_, 

I 

then  the  recombination  rate  when  n^  =  n_  ®  n  will  be 

2  2 

n  (s^  s„  -  84.8-) 


Hencff 


a  “  +  8^  -  s^s^> 


If  we  write  s^  and 's.  as 

=^=.l..^'o/V 

(5)  s_  -  1  - 


^  L_  =  L  mean  free  path  length 


TJien  for  high  pressures 


8^,  +  s_  -  s^8_  =  1 

_  2 
a  a  VTJTq 


and  for  low  pressures 


s .  +  8.  -  8. 8_  ^  —£ 


Consider  the  value  of  r. 


(8)  r^  =  4.05  X  10"^  (^)  om 

Hence,  for  the  large  nuclei  (R  Jk  r-)  the  value  of  the  recombination  coefficient 


given  by 


-  2 
P  a  VT® 


is  suggested.  It  will  be  assumed  that  equation  (9)  gives  the  small  ion  -  large 
nuclei  recombination  coefficient. 

For  nuclei  with  radii  R  ^r^  a  combination  of  the  Thomson  coefficient  and 
the  above  is  necessary.  For  equilibrium  conditions,  equation  (1)  then  becomes 


■[;5!  ' 


2PN  = 


if  en  is  small. 


where  n^  is  the  equilibrium  number  density  if  no  Altken  nuclei  are  present. 

Thus,  the  term  iPN  is  important  when  any  deficit  in  the  small  ion  density 


“6i>«  c 

ocoiirs.  Lik^Hlae,  any  jxeesa  obaeryi^  in  the  small  ion  density  is  to  be  attri- 

bated  to  sources  other  than  cosmic  radiation.  The  possibility  bf- variations--,  - 
.  „  ■  .  "J  f-  :  ■:  .  "  ,1  .. 

ooboring  in  the  cosmic  Inay  prodnctlon.:  are  remote  at  the  jj^titudes  considered 

except  on  special  occasions  (solar  flares  etc.).  Since  the  advent  of  naolear 

<1 

testing  and  in  viev  of  the  residence  time  of  radioa(^tive  material  near  and  . 
above  the  tropopause  (List,  1962),  an  excess  of  small  ions  i,8  expected  and  is 
observed  at  the  base  of  the  stratosphere. 

Consider  the  data  of  flight  58?,  figure  l6,  tabulated  in  Table  IV.  If 
one  divides  the  termISN  by  the  mean  .ionic  velocity  one  obtains  the  geometrical 
cross  sectional  area  of  the  nuclei  per  cm^  if  they  have  radii  R^rrQ.  If  the 
nuclei  in  the  main  have  radii  R«  r^,  the  area  obtained  is  fictitious.  It 
should  be  noticed  though  that  as  far  as  recombination  is  concerned  the  larger 
nuclei  are  weighted  according  to  their  cross  sectional  area.  Moreover,  the 
large  nuclei  (R  s.  r^)  become  increasingly  more  important  at  lower  pressures 
since  p  is  not  pressure  dependent  whereas  a  is  proportional  to  the  presshre. 

The  nuclei  (R<r^)  are  associated  with  a  recombination  coefficient  which  is 
not  nearly  as  pressure  dependent  as  a  save  for  very  small  R.  The  'excess* 
entries  in  Table  IV  under  2  PN  arise  since  here  n>njj.  In  this  region.  It  must 
be  assumed  that  cosmic  rays  are  not  the  only  important  production  factor.  The 
quantity  of  radioactive  material  available  for  small  ion  production  has  been 
measured  by  Machta  and  List  (1959)  and  List  (1962), 

One  can  at  best  say  that  this  is  a  variable  quantity.  But  to  demonstrate 
that  it  Is  present  in  sufficient  quantity  to  affect  the  ion  density  note  figure 
17.  On  this  flight  there  is  a  pronounced  layer  with  high  ton  density  restricted 
to  the  pressure  range  150  mb  -  110  mb.  This  layer  is  probably  in  adveotive 
motion  descending  through  the  stratosphere.  The  descending  motion  is  suggested 
by  the  temperature  profile  associated  with  the  ion  peak.  The  ion  density  in 
this  layer  approaches  twice  the  level  outside  the  layer.  An  additional  produo- 
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TABLE  IV  ji 


p(mb) 

n(3£S£) 

3 

cm 

-  3  i _ : 

‘  0  (saC) 

-  ap® 

sec 

800 

850 

0 

CO 

0^ 

2.1x10“^ 

.58x10“^ 

500 

2500 

2000 

2.1 

excess 

300 

4300 

3220  . 

excess 

150 

5800 

4160 

1-9 

excess 

70 

5000 

4390 

1.1  * 

excess 

30 

3750 

4120 

.53 

.46ri0“^ 

10 

^  2500 

3780 

,18 

.90x10"^ 

tion  rate  which  la  approxiinately  twice  that  of  cosmic  rays  must  be  assumed. 
This  layer  is  analogous  to  the  ozone  rivers  discussed  in  the  previous  section 
on  ozone. 

The  lower  IJjait  of  the  particle  size  responsible  for  the  deficit  of  small 

ions  at  high  altitude  is  partly  determined  by  the  collector  used  in  the 

measurements.  Since  the  collector  will  saturate  on  lonS  with  a  mobility  of 
2  • 

1  22! — at  H  T  P,  it  will  saturate  w.r.t.  ions  of  about  l/lOO  of  this 
volt  sec 

mobility  at  10  mb.  The  mobility  of  large  ions  at  the  surface  is  about 

5  X  10**^  °i7  -  ,  These  have  a  radius  of  about  5  x  10“^  cm  end  will  not  be 
volt  sec 

collected  at  high  altitude.  The  intermediate  ions  have  a  mobility  of  about 

7  X  10“^  — ^  —  at  N  T  P  and  a  radius  of  about  4  x  10“"^  cm.  These  will  be 
volt  sec 

collected  at  high  altitudes.  Hence,  it  appears  that  the  particles  must  be  on 
the  order  of  10“^  cm  radius  or  larger. 

To  get  an  idea  of  the  magnitude  of  the  nuclei  number  density  required  to 
explain  the  observations,  consider  the  deficit  pi*esent  at  10  mb  in  flight  587 
given  in  Table  I.  Inserting  a  mean  velocity  of  the  small  ions  of  3  x  10^  -|S 


ve  get 


=  3  X  :10  csr 


If  all  of  the  particles  have  a  radixis  ofi|  5.2  x  lo“®  Cm,  the  nainber  density 
required  is  about  350  per  cm'^.  If  the  radius  is  inoreasiSd  by  a  factor  of  lO, 
the  number  density  needed  is  only  3*5  per  cm^.  Hence,  contrary  to  ■jfhe  state¬ 
ment  of  Junge,  Chagnon  and  Manson  (1961),  the  ion  density  measurements  are  more 

than  sufficiently  sensitive  to  particulate  natter.  The  concentrations  of  dust 

1  /■ 

given  by  the  ion  density  measurements  are  at  least  a  factor  of  20  higher  than 
those  given  by  Junge.  In  view  of  the  uncertainties  and  inherent  difficulties 
in  the  technique  used  by  J\xnge,  however,  it  is  concluded  that  the  ion  density 
measurements  are  more  representative  of  ambient  conditions.  It  should  be  men¬ 
tioned  that  the  data  point  chosen  for  the  calculation  and  comparison  is 
typical. 

The  data  shown  in  figures  16  and  18  are  nearly  identical  to  the  profiles 
obtained  for  about  of  the  flights  made  since  this  program  was  initiated  in 
1958.  However,  there  ai*e  oases  where  the  nuclei  content  of  the  atmosphere  is 
very  high  (figure  17).  Before  proceeding,  the  following  question  must  be 

I 

answered.  If  nuclei  are  so  effective  for  removing  small  ions  at  high  altitude 
where  their  number  density  is  so  low,  then  why  are  there  any  small  ions  in  the 
troposphere  where  the  nuclei  count  is  high?  The  answer  to  this  question  is 
the  following I  At  high  pressures  (p^l50  mb)  the  Thomson  recombination  coef¬ 
ficient  (o)  is  as  large  as  ^  for  nuclei  with  Hence,  for  the  Aitken 

nuclei  term  to  be  Important,  the  Aitken  nuclei  must  exist  in  quantities  com¬ 
parable  to  the  small  ions.  This  is  not  the  case  since  N  is  only  about  200 
per  em^  in  the  troposphere.  Furthermore,  the  larger  nuclei  (R>>rQ)  are  lost 
rapidly  from  the  troposphere  by  washout,  whereas  the  small  nuclei  (R  eg  r.) 


require  larger  supersaturation  in  order  for  condensation  to  occur  on  them.  On 
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ij  .  ifi!  ’i 

ooo^lon  ttipr  should  be  prefent 


s^f^eient  quantit^i  however,  and  thlq  is 
“boiTi  out  by  figure  1?  and  <&ta  presented  elsewhere  (Kroe^ngt  I960),  Moreover,, 

^  if'  ■  .  .  .  '  ,  '  ' 

at  very  low  altitude  in  the  troposphere,  natural  radioactive  material  is  • 

"  ’  ^  s',  “  I  t  ,,  ' 

throwp  up^ilohg  vith  largb  nuclei,  thus  masking  their  presence. 

The  |erm5pil/v  has  another  interesting  aspect.  Since  this  is  the  geo¬ 
metrical  cross  section  of  the  particles .  per  cm^  if  ,iR  fe  r^,  one  might  e:^eot, 

/  i 

on  this  basljSu  a  variation  in  the  atmospheric  transparency.  Zacharov  (1952) 
indeed  has  found  such  a  variation  by  analjr^ing  the  10-3rear  measurements  of 
atmospheric  transparency  by  Abbot  and  Fowls  (1913)  and  has  attributed  it  to 
particulate  matter  of  a  similar  size  and  quantity  suggested  by  the  ion  density 
measurements. 

For  example,  if  the  nuclei  are  distributed  more  or  less  unifonnly  through¬ 
out  an  atmospheric  scale  height  (10  km),  we  would  estimate  a  total  geometrical 
cross  sectional  area  of  ,03  om^ .  Thus,  an  attenuation  of  35^  in  the  vertical 
would  result.  This  assumes  that  the  attenuation  cross  section  is  the  same  as 
the  geometrical  cross  section.  This  is  so  for  particles  near  the  wavelength 
of  light.  Zacharov’ s  attenuation  factor  was  wavelength  dependent  as  indicated 
belowt 
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tdiere  t  is  the  transparency  of  the  atmosphere.  From  this  selective  attenuation, 
Zacharov  concluded  that  the  particles  could  not  have  diameters  greater  than 
10  cm.  Giovanelll  (1954),  on  the  other  hand,  using  Zacharov' s  data  concludes 
that  the  particles  could  have  diameters  larger  than  10"^  cm.  Qlovanelli 
estimates  the  total  number  of  particles  to  be 
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He  farther  assime^  t^t  they  are  dl8tribut«>d  throughout  a  10 -km  range 
and  arrives  at  a  nuniber  density  of  3  particles  jber  om^.  This  is»  for  pradti- 
oal  "purposes,  the  sane  as  that  required  by  thel  ion  density  measurements J  '  It 
is  eoncluded  that  the  data  of  Junge,  Chagrion  and  Hanson  (I961)  are  in  error  ai 
these  high  altitudes.  There  are  several  oritioal  points  to  be  made  regarding 
thd  data  of  Junge,  Chagnon  and  Hanson. 

1.  At  the  operating  temperature  of  the  expansion  lohamber  (20®e) 
used  to  determine  the  nudlei  count,  the  vapor  pressure  of  water  is  23*4  sib.  ’ 
This  pressure  corresponds  to  an  altitude  of  25  hm  in  the  atmosphere.  ITrom 
Junge *8  description  of  the  apparatus,  it  is  unclear  as  to  how  the  measurements 
could  be  meaningful  above  25  km.  Prior  to  each  expansion  of  the  ohasiber 

1000  00  of  ambient  air  were  supposedly  drawn  through  the  system  within  2tJ 
seconds  via  an  inlet  and  outlet  valve.  A  pump  was  connected  to  the  outlet. 
However,  an  experiment  was  performed  by  the  author  to  determine  the  evaporation 
rate  of  water  at  this  oritioal  pressure  (23.4  mb).  In  25  seconds  a  1  square 
inch  wetted  piece  of  black  velvet  at  20®C  would  evaporate  3500  oe  of  water 
vapor.  From  this  one  can  only  deduce  that  no  fresh  air  entered  the  chamber  and 
only  that  water  vapor  left  both  inlet  and  outlet  valves  the  monent  they  were 
opened  for  sampling. 

2.  The  efficiency  of  the  detectors  was  calculated  and  not  determined 
^erimentally.  This  is  Just  as  serious  as  the  previous  criticism  except  that 
it  applies  at  all.  altitudes. 

3.  No  descent  data  were  obtained  since  the  equipment  always  failed 
at  or  before  reaching  ceiling  altitude, 

4.  The  statistics  involved  in  some  of  the  measurements  was  poor. 

How  in  regard  to  the  ion  density  measurements  there  is  only  one  known 

factor  which  could  affect  the  results.  That  is  the  electrostatic  effect  of  the 


n 


balloon.  The  magnitude  of  this  effect  is  shown  In  f 16  where  both  asjjsent 
and  descent  data  are  presented  for  t^ie  sane  fli'lht.'  itie  variation  bet^ah  ascent 
and  descent  at  10  mb  is  a  considerable  factor  and  it  is  on^of  /'the  largest  ob¬ 
served.  Moreover)  it  does  not  affect  the  above,  calculatloiis  sjensibly.  In  all 
of  the  soundings  the  package  is  suspended  200  ft.  below  the  b|lloon  to  i^nlnlze 
this  effect. 

Exception  to  the  above  stataaent  is  allowed  for  the  two> flights  shown  in 
figures  20  and  21.  On  these  two  flights  and  only  on  these  two*  there  was  the 
possibility  that  an  electric  field  existed  between  the  outside  of  the  collector 
and  the  gondola .  This  external  electric  field  would  have  the  effect  of  reducing 
the  collection  rate  of  negative  ions.  These  two  flights  are  presented  mainly 
for  this  purpose.  Yet,  there  is  the  possibility  that  they  do  represent  ambient 
conditions  since  the  flights  followed  the  Perseld  and  Arletid  meteor  showers . 
Furthermore,  the  sporadic  meteor  Influx  had  already  reached  its  maximum  prior  to 
these  two  flights. 


Fig.  20,  The  ion  density  and  teinperature  nrofile  of  flight  602  for  Jnly  13»  1962, over  Minneapolis. 
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Fig.  21.  The  ion  density  and  temperature  profile  of  flight  603  for  July  24,  1962,  over  Minneapolis. 


Fig.  21.  Schematic  diagram  of  the  electrometer  used  in  the  ion  density  measurements. 
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